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Abstract 

In this paper we establish the existence of a square integrable occupation 
density for two classes of stochastic processes. First we consider a Gaussian 
process with an absolutely continuous random drift, and secondly we handle the 
case of a (Skorohod) integral with respect to the fractional Brownian motion with 
Hurst parameter H > ^. The proof of these results uses a general criterion for 
the existence of a square integrable local time, which is based on the techniques 
of Malliavin calculus. 
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1 Introduction 



Local times for semimartingales have been widely studied. See for example the mono- 
graph [11] and the references therein. On the other hand, local times of Gaussian 
processes have also been the object of a rich probabilistic literature; see for example 
the recent paper by Marcus and Rosen. A general criterion for the existence of a 
local time for a wide class of anticipating processes, which are not semimartingales 
or Gaussian processes, was established by Imkeller and Nualart in [6]. The proof of 
this result combines the techniques of Malliavin calculus with the criterion given by 
Geman and Horowitz in This criterion was applied in [6j to the Brownian motion 
with an anticipating drift, and to indefinite Skorohod integral processes. 

The aim of this paper is to establish the existence of the occupations densities 
for two classes of stochastic processes related to the fractional Brownian motion, using 
the approach introduced in |6]. First we consider a Gaussian process B = {Bt,t G 
[0, 1]} with an absolutely continuous random drift 



where u is a stochastic process measurable with respect to the cr-field generated by 
B. We assume that the variance of the increment of the Gaussian process B on an 
interval [s,t] behaves as \t — s\'^^, for some p € (0,1). This includes, for instance, 
the bifr actional Brownian motion with parameters H,K G (0, 1). Under reasonable 
regularity hypotheses imposed to the process u we prove the existence of a square 
integrable occupation density with respect to the Lebesque measure for the process 



Our second example is represented by the indefinite divergence (Skorohod) 
integral X = {Xt,t G [(0,1]} with respect to the fractional Brownian motion with 
Hurst parameter H G (^,1), that is 



We provide integrability conditions on the integrand u and its iterated derivatives in 
the sense of Malliavin calculus in order to deduce the existence of a square integrable 
occupation densities for X. 

We organized our paper as follows. Section 2 contains some preliminaries on 
the Malliavin calculus with respect to Gaussian processes. In Section 3 we prove the 
existence of the occupation densities for perturbed Gaussian processes and in Section 
4 we treat the case of indefinite divergence integral processes with respect to the 
fractional Brownian motion. 




X. 




2 



2 Preliminaries 



Let {Bt,t € [0, 1]} be a centered Gaussian process with covariance function 

R{t,s) :=E{BtBs), 

defined in a complete probability space {0,,J^,P). By Ti. we denote the canonical 
Hilbert space associated to B defined as the closure of the linear space generated by 
the indicator functions {l[o^j],t G [0, 1]} with respect to the inner product 

(l[0,t],l[0,s])w = R{t,s), s,t G [0,1]. 

The mapping 1[q^^] Xf can be extended to an isometry between 7i and the first 
Gaussian chaos generated by B. We denote by B{ip) the image of an element ip £ 7i 
by this isometry. 

We will first introduce some elements of the Malliavin calculus associated with 
B. We refer to [10] for a detailed account of these notions. For a smooth random 
variable F = f {B (ipi) , . . . , B (ipn)) , with (/jj G H and / G C^iR"") (/ and all its 
partial derivatives are bounded) the derivative of F with respect to B is defined by 




For any integer k > 1 and any real number p > 1 we denote by D'^'^ the Sobolev space 
defined as the the closure of the space of smooth random variables with respect to 
the norm 

i=i 

Similarly, for a given Hilbert space V we can define Sobolev spaces of ^-valued random 
variables B'''P{V). 

Consider the adjoint 6 oi D in L^. Its domain is the class of elements u G 
L^{n;n) such that 

E{{DF,u)n)<C\\F\\2, 
for any F G B^'^, and 5 (n) is the element of L'^{Q) given by 

E{6{u)F) = E{{DF,u)n) 

for any F G B^'^. We will make use of the notation 6{u) = Jq UsSBg. It is well-known 
that ro^'^(7^) is included in the domain of 5. Note that E{6{u)) = and the variance 
of d{u) is given by 

E{d{uf) = E{\\u\\l^) + E{{Du, (Du)* 
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if n E ©^'^(W), where (Du)* is the adjoint of Du in the Hilbert space Ti (^7i. We 
have Meyer's inequahty 

for any p > 1. We will make use of the property 

F5{u) = 5{Fu) + {DF, u)n- (2.3) 

if F € B^'^ and u G Dom((5) such that Fu G Dom((5). We also need the commutativity 
relationship between D and 5 

D5{u) = u+ [ Dus6Bs, (2.4) 

if n G B-'^'^(7^) and the process {DgU, s G [0, 1]} belongs to the domain of 5. 

Throughout this paper we will assume that the centered Gaussian process 
B = {Bt,te [0,1]} satisfies 

Ciit - sfP < E{\Bt - B,\^) < C2{t - sfP, (2.5) 

for some p G (0, 1) with Ci,C2 two positive constants not depending on t,s. It will 
follow from the Kolmogorov criterium that B admits a Holder continuous version of 
order 5 for any 6 < p. 

Throughout this paper we will denote by C a generic constant that may be 
different from line to line. 



Example 1 The bifractional Brownian motion (see, for instance J^), denoted by 
B^'^ , is defined as a centered Gaussian process starting from zero with covariance 

Rit,s) = ^{{t'^ + s^^f-\t-sr^) (2.6) 

where H G (0, 1) and K G (0, 1]. When K = 1, then we have a standard fractional 
Brownian motion denoted by B^ . It has been proven in that for all s < t, 



<2i-^|t-s|2^^^ (2.7) 



so relation \2. 5|] holds with p = HK . A stochastic analysis for this process can be 
found in J7] and a study of its occupation densities has been done in f^, flEI- 
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For a measurable function x : [0, 1] — > M we define the occupation measure 

fj.{x){C) = / lcixs)ds, 
Jo 

where C is a Borel subset of R and we wiU say that x has on occupation density with 
respect to the Lebesque measure A if the measure /i is absolutely continuous with 
respect to A. The occupation density of the function x will be the derivative For 
a continuous process {Xt,t € [0, 1]} we will say that X has an occupation density on 
[0, 1] if for almost all u; G fi, X(uj) has an occupation density on [0, 1]. 

We will use the following criterium for the existence of occupation densities 
(see [6]). Set T = {{s, t) G [0, 1]^ : s < t}. 

Theorem 1 Let {Xt,t € [0,1]} be a continuous stochastic process such that Xt G 
p2,2 every t G [0, 1] . Suppose that there exists a sequence of random variables 
{Fn,n > 1} with Uni-^" 7^ 0} = r2 a.s. and Fn G B^'^ for every n> 1, two sequences 
an > 0,6n > 0, a measurable bounded function 7 : [0, 1] M, and a constant 9 > 0, 
such that: 

a) For every n > 1, \t — s\ < 5n, and on {Fn 7^ 0} we have 

{jD{Xt-Xs),l[s,t])H> an\t-sf, a.s.. (2.8) 

b) For every n > 1 

[ E{{^DFn,l[s,t])H)\t - s\-^dtds < 00. (2.9) 
Jt 

c) For every n > 1 

^\t-s\~'^''dsdt< 00. (2.10) 

Then the process {Xt,t G [0,1]} admits a square integrable occupation density on 
[0,1]. 

Remark 1 The original result has been stated in JFJ with 9 = 1 in the case of the 
standard Brownian motion. On the other hand, by applying Proposition 2.3 and 
Theorem 2.1 in it follows easily that this criterium can be stated for any 9 > 0. 



/ E(^Fn{^'^''DD{Xt-Xs),lfs% 
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3 Occupation density for Gaussian processes with ran- 
dom drift 



We study in this part the existence of the occupation density for Gaussian processes 
perturbed by a absolute continuous random drift. The main result of this section is 
the following. 

Theorem 2 Let {Bt,t € [0, 1]} be a Gaussian process satisfying ^2. 5|) . Consider the 
process {Xt,t G [0, 1]} given by 

Xt = Bt+ / Usds, 
Jo 

and suppose that the process u satisfies the following conditions: 
1. nGB2'2(L2Qo,l])). 

^- ^ ((/o^ 11-^^^*11 '^^y'^^ < oo, for some q>l,p> jh^. 
Then, the process X has a square integrable occupation density on the interval [0, 1]. 



Proof: We are going to apply Theorem[TJ Notice first that Xt E B^.s fQj^ g,!! t G [0, 1]. 
For any < s < t < 1, using (j2.4p and (|2.5p we have 



s,t] 



H 



> Ci{t-s 



,2p 



Dur-dr, 



H 



> Ci{t - sfP - y^{t - sY \\Dur\\-^^dr. 



By Holder's inequality, if | + ^ = 1, we obtain 



DurWf^ dr < (t — s) 



iDurWj^ dr 



Fix a natural number n > 2, and choose a function ipn{x), which is infinitely differ- 
entiable with compact support, such that <^n(x) = 1 if |x| < n — 1, and (fnix) = 0, if 
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\x\ > n. Set Fn = ^Pn Jq \\Dut\\^ dtj ^ j . The random variable F„ belogs to D^''^. 
In fact, it suffices to write -F„ = (pniG), where 

[Dur, hf)-^ dr, 



G = sup 

heL9([0,i];W) -'0 
\M<1 



which imphes 



yjG)DG\\^ < 



^ II ' I 



sup 

/ieL''([0,l];H) 

1 

V 



n 



Then, on the set {F„ 7^ 0}, ^/J ||DMt||^ dtj ^ < n, and we get 

= {t - s)^P \Ci - ny^{t - s)-^-^ 

and property a) of Theorem [1] holds with a suitable choice of an and 5n because 
- — p > 0, and with 6 = 2p and 7 = 1. 



Finally, conditions b) and c) can also be checked: 



E 



[DFn,l[s,t]) 



H 



|t-s|2P 



dsdt < \fC2 



E{\DFn\\n) 
\t - s\P 



dsdt < 00, 



and 



E 



Fn{D\Xt-Xs),l' 



s,t] 



dsdt 



< \\Fn\LC2 



\t-s\^P 

E{\\D^{Xt-Xs)y^,) 

— dsdt < 00, 



because 



E{\\D\Xt-Xs)\\^^,) = I'd 

< {t-s)^E 



'Urdr 



< / E(\\D^Ur\\..^o)dr 



I F) Ur 1 1 •^02 dv 
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and | — 2p=l — i — 2/9> —1, because p > 2(i-p) ■ 



Remark 2 These conditions are intrinsic and they do not depend on the structure of 
the Hilbert space 7i. In the case of the Brownian motion, this result is slightly weaker 
than Theorem 3.1 in 16^, because we require a little more integrability. 



4 Occupation density for Skorohod integrals with re- 
spect to the fractional Brownian motion 

We study here the existence of occupation densities for indefinite divergence integrals 
with respect to the fractional Brownian motion. Consider a process of the form Xt = 
JqUsSB^^, t £ [0,1], where B is fractional Brownian motion with Hurst parameter 
H G (i,l), and u is an element of I])^''^{L'^{[0,1])) C Dom((^). 

We know that the covariance of the fractional Brownian motion can be written 
as ^ 

E{B^Bf)= [ r ct>{a,(5)dadf], (4.1) 
Jo Jo 

where (/>(«, /3) = H{2H - l)\a - For any < s < t < 1, and a G [0, 1] we set 

/,,t(a):= [ (t>{a,(3)dadp. (4.2) 

J s 

We also know (see e.g. |10j ) that the canonical Hilbert space associated to B satisfies: 

l2([0,1])cL^([0,1])cW. (4.3) 

The following is the main result of this section. 

Theorem 3 Consider the stochastic process Xt = /q Us5B^ where the integrand u 
satisfy the following conditions for some q > y~h o-ndp > 1 such that |+2 < H{p+1): 

11) u G]D)3'2(l2([0,1])). 

12) /J f^[E{\DtUs\n + ^(11 \DtDus\ 11^) + ^(11 \DtDDus\ ||?:,^^)]dsdt < oo. 

13) E ({utr^^'^''^^^^ dt < oo. 

Then the process {Xt,t G [0,1]} admits a square integrable occupation density on 
[0,1]. 
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Proof: We will use the criteria given in [S] and recalled in Theorem [TJ Condition 
II) implies that Xt G B'^'^ for all t € [0, 1]. On the other hand, from Theorem 7.8 in 
[7] (or also by a slightly modification of Theorem 5 in [T]) we obtain the continuity 
of the paths of the process X. Note that from Lemma 2.2 in [6] corroborated with 
hypothesis 13). we obtain the existence of a function 7 : [0, 1] {— Ij 1} such that 
'jtUt = \ut\ for almost all t and uj. 

We are going to show conditions a) , b) and c) of Theorem [H 

Proof of condition a): Fix < s < t < 1. From (j2.4[) we obtain 

D{Xt - Xs) = ^/l[,,t] + / Dur6B^, 

J s 

and we can write 

{^{Xt-X,), = (|n|l[,,t],lM)w + (7 f Dur6B^,l[,^t])n. (4.4) 

J s 

We first study the term 

{\u\l[s,t],Ms,t])H = / \ua\(l){a, f3)dadp = \ua\fs,ti(^)da- 



For any g > 1 we have 

£;(|5f-i?f|2) = jj,^t{a)da 

q _ q 

( ha I («))«+' (|ua|/s,t(a)) fs,t{(y)da, 

and using Holder's inequality with orders and g + 1, we obtain 

t g J 1 

E{\B^ -B^\^) < (^^ |n,|/,,t(a)da^ (^^ |u„rVs,t(a)da^ . 

Hence, using that 

/s,t(a) < /o,i(«) = H{2H - 1) \a - (3\^''~Hp = H (a^^-^ + (1 - af""-') < H, 

Jo 

we get 

2H(q+l) _ 1 

\ua\fsA<^)da > C\t - s\^^ Zg " , (4.5) 
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where Zq = /J^ \uc\~'^da. 

On the other hand, for the second summand in the right-hand side of (j4.4p 
we can write, using Holder's inequahty. 



DaUrSB^ 



fs,t{a)da 



< 



£-1 
P 



DaUrSBj. 



da 



We can write 



fs,t{a)p~^da 



p-i 

p 



= CH 
< CH\\l[s,t]* 



I i2H-2 



(4.6) 



(4.7) 



where ch = H(2H — 1). Young's inequahty with exponents a and b in (l,oo) such 
that 1 + 1 = 2 - i yields 



I Ifs /-I * I ■ 1^"^ nil ^ l|l 



s,t]\\L''( 



\2H-2 



1[-1,1]| 



(4.8) 



Choosing b < and letting t] = ^ < 2H - ^ we obtain from (l46l) . (fiTfl) . and (|M 



7 / Dur6B^,l[st] 

J s 



n 



< C\t-s\ 



DaUrdBr 



da 



Now we will apply Garsia-Rodemich-Ramsey's lemma (see fSj) with $(x) = x^, ^(a;) : 

m+2 

j; p and to the continuous function Us = Jq DaUrdBr (use again Theorem 5 in [1] 
and we get 



DUr6Br,'yl[s^t] 



n 



<c\t-sr^Y^,p, 



(4.9) 



where 



Yr 



"1 rl A I fy 



m,p 



Jo Jo 



f^DglirdBrf 



dxdyda. 
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Substituting and into ([O]) yields 



2g(9 + l) -i mi 



2g(9+l) 



1 1 



\t — s\ 1 I Zq ' — C|t — s\ Ym,p I , 



where 5 = V ~^ ^ ~ — With a right choice of rj the exponent 5 is positive, 
provided that — | — ^ > 0, because rj < 2H — |. Taking into account that 
^ < 1 - H, it suffices that 

m> - + 1- H. (4.10) 

P 

We construct now the sequence {Fn,n > 1}. Fix a natural number n > 2, 
and choose a function ipn{x), which is infinitely differentiable with compact support, 
such that (fnix) = 1 if < n — 1, and (pn{x) = 0, if > n. Set F„ = {G ), 
where G = Zg + Y^^p- Then clearly the sequences a„ and 5n required in Theorem [1] 
can be constructed on the set {Fn / 0}, with = 2H + 

It only remains to show that the random variables Fn are in the space B^'^. 
For this we have to show that the random variables HDZ^Il^ and |j 1)1^,^11-;^ are 
integrable on the set {G < n}. First notice that, as in the proof of Proposition 
4.1 of [6], we can show that E {\\DZq\\^) < oo. This follows from the integr ability 
conditions 13) and 

/ E {\\Dutf^) dt < oo, (4.11) 
Jo 

which holds because of 12), the continuous embedding of of -L"h([0, 1]) into Tl (see 
^J), and the fact that pH > 1. On the other hand, we can write 

Jo Jo Jo 
where £,x,y,a = Jy DaUr5Br- Thus 

\\F>Ym,p\\H < P / / \^x,y,a\'^~^ \\D^:r,y,a\\n\x - y\~"'~'^dxdyda 
Jo Jo Jo 

< p{Ym,p)^ {^j^ ^ ^ \\Dix,yAH\^-y\~""~^dxdydo^ . 
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Now, to show that l(G<n) ||-Dl^,p||>^ belongs to L^{Q), it suffices to show that the 
random variable 



Y = 



C C [' m.,y,a\fn\^-y\-"'-^dxdyda 
Jo Jo Jo 



has a finite expectation. Since, for any < y < x < 1 
we have 

Y < C 



m\\l[y^r,]Dau\\'^\x - y\ "^dxdyda 
+ ! I I \\ ! DDaUs6B^\\^\x-y\-'^'^dxdyda) 

Jo Jo Jo Jy J 

C{Y^ + Y2). 

From the continuous embedding of ^^([0, 1]) into TL, we obtain 



1 fi fi 




Yi < C 1^ \\liy^^]Dauf^^^H(^[Q^^)\x - y\ ™ ^dxdyda 



< C\x-y\P^^^ f C C r \D^Urf\x-y\-'^-'^drdxdyda 

Jo Jo Jo Jy 



Hence, E(Yi) < oo, by Fubini's theorem, Proposition 3.1 in [6j and condition 12), 
provided 

m<pH-l. (4.12) 
On the other hand, using the estimate (12. 2p . and again the continuous embedding of 
L"^([0, 1]) into n yields 



DD^Us6B^\\p^ < C^(||Z),i?n.l[j^,,](.)||^«2 + ||Z)aZ?i?u.l[j,,:,](-)| 

< C \\\Dc,Du.\ l[y,x](-)||ii///([o,i];^^) 

+ III-Dq-D-Du.! l[y,a.](-)||^i/H([o,l];'H®2) ) 

< Clx-ylP""-'!^!^^ E{\\\D^Dur\r^)dr 

-f 



IP 



+ ^ E [\\\Da,DDur\\\l^^2) dr 
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As before we obtain £(¥2) < 00 by Fubini's theorem and condition 12), provided 
(j4.12p holds. Notice that condition ^ + 2 < H{p + 1) impHes that we can choose 
an m such that and ([il^ hold. 

Proof of condition b): Define = {G < n}. Then, condition b) in Theorem [T] 
follows from 

/ E{{jDFn,lis,t])H)\t-s\-'dtds < C [ E{lA„\{lDG,l[,^t])n\)\t-s\-'dtds 

JT JT 

< C E {1a jDGWn) [ \t-s\^-^dsdt < 00, 

JT 

since E WDGWn) <oo8.ndd-H = H + ^<l. 
Proof of condition c): We have 



Hence 



Jl{s,t) + J2{s,t) + J3{s,t). 



For i = 1, 2, 3, we set 

Ai = E {En jjt - s\-'^^ \Ji{s,t)\dsdt^ . 



Let us compute first 



Ai < G \t-s\^''-^'E{\\\D^up\ly,^t^{(3)\\H^^)dsdt 

JT JT 

^ ^ ItIt^^~ '''''''' y^'^^'^y) ' 



where 

"1 fi 



ip{l3,y)= / E {\DaUi3\\D:,Uy\) (j){a,x)(j){P,y)dadx. 

JO JO 
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By Fubini's theorem Ai < oo, because 2H — 29 > —2, which is equivalent to q > H, 
and 



1 fi 



^0 



^{P,y)dpdy<E(\\\Du\\\l^^2 



and this is finite because of the inclusion of L'^{[0, 1]) mTi (|4.3p . In the same way we 
can show that ^2 < oo- Finally, 



A. 



E R 



t-s 



-29 



T JT 



7 



J s ' 



dsdt 



< C 



\t-s\''H-^'E 



T JT 



DDur5B^ 



dsdt, 



and we conclude as before by using for example the bound (12. 2p for the norm of the 
Skorohod integral and the condition 12). ■ 



Remark 3 If p 



1+V17 



, then i + 2 < H{p + 1) for all H > \. 
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